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The hrj /Et signal at the ongoing LHC experiments is simulated with Pythia in 
the mSUGRA and other models of SUSY breaking. Special attention is given on the 
compatibility of this signature with the low mass neutralino dark matter (LMNDM) 



^ i' scenario consistent with WMAP data. In the mSUGRA model the above signal as 

^ '. 
,^ , well as the LMNDM scenario are strongly disfavored due to the constraints from the 



on going SUSY searches at the LHC. This tension, however, originates from the model 
dependent correlations among the parameters in the strong and electroweak sectors 
of mSUGRA. That there is no serious conflict between the LMNDM scenario and the 
LHC data is demonstrated by constructing generic phenomenological models such that 
^^ \ the strong sector is unconstrained or mildly constrained by the existing LHC data and 

parameters in the electroweak sector, unrelated to the strong sector, yield DM relic 
density consistent with the WMAP data. The proposed models, fairly insensitive to 
j_j ■ the conventional SUSY searches in the jets + ^t and other channels, yield observable 

signal in the suggested channel for C ~ Ifb of data. They are also consistent with 
the LMNDM scenario and can be tested by the direct dark matter search experiments 
in the near future. Some of these models can be realized by non-universal scalar and 
gaugino masses at the GUT scale. 
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1 Introduction 

Proton - proton collisions at the LHC are now producing data at a center of mass energy 
(y/s) = 7 TeV. There is no evidence of any new physics beyond the standard model (SM). 
However, it has already been shown by the ATLAS [1] and CMS [2] collaborations that even 
with a small integrated luminosity (£) = 35 pb~^, supersymmetry (SUSY) can be probed 
much beyond the existing limits on the sparticle masses obtained by the LEP [3] or Tevatron 
[4] experiments. 

The negative results of new particle searches at the LHC have been interpreted in terms of 
the simplest gravity mediated SUSY breaking model - the minimal supergravity (mSUGRA) 
[5] model- which has only five free parameters including soft SUSY breaking terms. These are 
rrio (the common scalar mass), mi/2 (the common gaugino mass), Aq (the common trilinear 
coupling parameter), all given at the gauge coupling unification scale (M^ ~ 2 x 10^^ GeV); 
the ratio of the Higgs vacuum expectation values at the electroweak scale namely tan/3 and 
the sign of fi. The magnitude of fi is determined by the radiative electroweak symmetry 
breaking (REWSB) condition [6]. The sparticle spectra and couplings at the electroweak 
scale are generated by renormalization group evolutions (RGE) of the above soft breaking 
masses and the coupling parameters. The non-observation of signal, in particular in the 
jets + /Er, channel leads to exclusion plots in the mo — mi/2 plane. A large number of 
phenomenological analyses have also addressed the issue of SUSY search at LHC-7 TeV 
experiments [7, 8]. 

In a hadron collider the dominant source of SUSY signals in the m / -|- n jets + /^t 
channel is the pair production of the strongly interacting sparticles - the squarks(g) and 
gluinos(^) - in various combinations. Throughout this paper / stands for e and /i unless 
stated otherwise. Thus the bounds from ATLAS [1] and CMS [2] primarily exclude some 
parameter space with relatively low rriQ and mi/2 which translates to certain combinations 
of squark q and gluino g masses ^. For example, the non-observation of the 0/ + jets + ^t 
signal implies that for nearly mass degenerate squarks and gluinos m~ ^ m^ > 775 GeV (see 
the second paper of [1]). Here m~ stands for the average mass of the L and R type squarks. 

On the other hand in several regions of the parameter space of the minimal supersymmet- 
ric extension of the SM (MSSM) with conserved R-parity but without specific assumptions 



"'The constraints become more severe due to the very recent £ = 1 fb ^ data as discussed briefly in 
Section 4. 



about the soft breaking parameters, the dark matter rehc density [9, 10] in the universe - 
low mass neutrahno dark matter (LMNDM) in particular - essentially depends on the prop- 
erties of the sparticles in the electroweak (EW) sector. This, e.g., is the case if the lightest 
neutralino, assumed to be lightest supersymmetric particle (LSP) {xi)y is bino like and all 
squarks are beyond the reach of the ongoing LHC experiments. 

However, due to the specific correlations among the sparticle masses in mSUGRA, the 
above bounds on m~ and m-^ would also impose stringent indirect mass bounds on the EW 
sparticles. This disfavors the LMNDM scenario. The bound quoted above from ATLAS 
data, e.g., implies mj ~ 398 and m~o ~ 125 ^. These model dependent bounds imply that 
the masses of the sparticles belonging to the EW sector are far above the direct lower limits 
from LEP [3] and too large for realising the LMNDM scenario. 

The observed dark matter (DM)relic density [Qh"^) in the universe [9, 10] has been 
precisely measured by the Wilkinson Microwave Anisotropy Probe (WMAP) collaboration 
and is bounded by 0.09 < Qh^ < 0.13 [11]. A possible mechanism of production of relic 
density in the above range in the LMNDM scenario is annihilation of a bino like LSP pair 
or bulk annihilation [10, 12, 13]. It may be recalled that, if all strongly interacting sparticles 
are heavy, relatively low mass neutralinos and R type sleptons (super partners of cr, hr 
and T/j) mainly contribute to this annihilation process. Coannihilation [13, 14] of a light 
neutralino and a nearly degenerate lighter stau mass eigenstate (ri) is another proposed 
mechanism for generating the observed relic density. The allowed LMNDM scenarios in the 
mSUGRA model, with emphasis on the above two processes, have been delineated in the 
figures in [15] using parameter spaces different from the conventional ones. In this paper we 
shall frequently refer to these figures. It seems that both the above processes are apparently 
in conflict with the recent LHC data. The incompatibility of DM relic density production 
by slepton coannihilation and the data from the LHC-7 TeV experiments have recently been 
noted in [16]. 

The indirect 'exclusion' of a light electroweak sector will have a bearing on direct detec- 
tion of DM [17] as well. The tension between the constraints obtained by the ongoing LHC 
experiments and the mSUGRA parameter space accessible to direct DM search experiments 
by the XENON [18] and the CDMS [19] collaborations have also been noted in the literature 



^Throughout this paper all masses, mass parameters and quantities having the dimension of mass are 
given in GcV unless stated otherwise. 



[16, 20]. Several groups have also reported on the prospect of constructing the mass of the 
neutralino by such experiments. It is estimated that if m~o ~ 150, then it might be possible 
to reconstruct this mass by measuring the energy spectrum of the recoiling nuclear targets 
[10, 21]. The recoil energy spectrum is insensitive to higher neutralino masses. Moreover the 
LMNDM scenario can be tested in an e"*" — e~ collider if the LSP mass is in the range 1-10 
[22]. However, from the results of direct DM search and/or various constraints from collider 
and astrophysical experiments it has been claimed that the above mass range is disfavored 
[23]. 

In view of the above discussions it is worthwhile to critically reexamine the constraints 
from the ongoing LHC experiments and their impact on LMNDM scenarios. This will be 
taken up in a later section. The main conclusion is that in view of the uncertainties in 
the data, some parameter space with low mo and mi/2 consistent with LMNDM cannot be 
conclusively ruled out. However, it must be admitted that there is a tension between the 
LHC data and the LMNDM scenario in the mSUGRA model. 

We remind the reader that before the advent of the bounds from the LHC, the bound 
rrifi > 114.4 on the lighter Higgs scalar mass {nih) from LEP [24] tightly constrained the 
parameter space with low niQ — mi/2, a part of which coincides with the parameter space 
corresponding to the LMNDM scenario. These constraints are particularly severe for low 
and intermediate values of tan/3 [25]. 

It was recently emphasized in [15, 26] that in order to revive the parameter space consis- 
tent with the LMNDM scenario a moderate to large negative values of the trilinear coupling 
(Aq) is called for. This is particularly important if tan/3 is not very large. In fact for Aq 
= and tan/3 = 3 - a choice frequently employed by the LHC and Tevatron experiments 
- the entire mSUGRA parameter space sensitive to the 35 pb~^ data is already excluded 
by the ruh bound. For sizable negative values of Aq the LMNDM scenario is realized for 
another reason. Here the lighter stau mass eigenstate (ri) becomes significantly lighter than 
the selectron or the smuon even for moderate tan/3 and LSP - fi coannihilation may occur 
efficiently. This happens for a minimum value of mi/2 much lower than the corresponding 
value for Aq = (see the figures depicting the parameter space allowed by WMAP data in 
[15]). 

This also gives rise to spectacular collider signatures as the lighter chargino (xi ) - or 
the second lightest neutralino (X2) dominantly decay into r rich final states. Predictions for 



experiments at LHC-14 TeV [15, 27] and LHC-7 TeV [8] were worked out. In this paper we 
shall restrict ourselves to non-zero values of Aq only. 

It is, however, well-known that the 0/ + jets + ^t signal is fairly insensitive to the choice 
of Aq and tan/3 (see the second paper in [1]). Throughout this paper we shall assume that 
the bounds obtained by the ATLAS and the CMS collaborations in this channel for fixed 
choices of these two parameters are valid for other choices as well. 

Another important consequence of non-zero Aq is that the lighter top squark mass eigen- 
state ti may be significantly lighter than the other squarks and could be copiously produced 
at the ongoing LHC experiments. They may come from two dominant sources; i) Direct 
ti — tl pair production and ii) g — )■ tit, if kinematically allowed. 

Our next task is to identify a signal which unlike the jets + /^t final state is sensitive 
to \Aq\. If tl is not the next lightest super particle, then the decay mode ti — > bxt may be 
its main decay channel resulting in final states rich in b-jets. If tan/3 is small (say, 5) then 
the electroweak gauginos will dominantly decay leptonically into e, /i or r channels with 
approximately equal probability. A viable signal in this case would be blj ^t [28] . 

For moderate or large tan/3 and non-zero Aq on the other hand the above EW gauginos 
decay dominantly into final states involving r's leading to very characteristic collider signals 
[8, 15, 27]. The price to be paid for r tagging efficiency may be adequately compensated by 
the large BRs of the EW gaugino decays into r rich final states. The main purpose of this 
paper is to study the observability brj ^t at the ongoing LHC experiments with emphasis 
on the LMNDM scenario. Here j is the number of jets in the signal and different choices of 
this variable will be considered. Occasionally, however, we shall also revisit the blj ^t signal 
[28]. 

As discussed above there is indeed a tension between the realisation LMNDM in 
mSUGRA and the preliminary data from the LHC. Should the experimental constraints 
be strengthened in future, the tension will further intensify. In view of this we propose 
a few generic models which are either mildly constrained or unconstrained by the current 
LHC data and are consistent with the low mass neutralino DM scenario. The main point 
is that the LHC data is sensitive to the masses of the strongly interacting sparticles while 
the realization of LMNDM in these models DM hinges on the properties of the electroweak 
sector. Thus if the two sectors are uncorrelated the above tension will cease to exist. These 
models are generic in the sense that their viability depends on certain mass hierarchies 



among the strongly interacting sparticles rather than on some specific choices of the masses. 
The important correlations among different mass hierarchies in SUSY models and the cor- 
responding collider signatures have been emphasized in the literature [29]. The parameters 
in the electroweak sector can be chosen independently. In fact all models where the sparti- 
cle masses in the EW sector consistent with the corresponding LEP limits, derived without 
assuming mSUGRA as the underlying model, are allowed in principle. These models are 
phenomenological in nature ^ although we shall comment on theoretical motivations wher- 
ever possible. Finally we shall discuss the possibihty of testing these models by the ongoing 
LHC experiments. 

The plan of the paper is as follows. In Section 2 we shall mainly concentrate on the 
phenomenological models unconstrained by the LHC data corresponding to £ = 35 pb~ 
above and assess the prospect of observing the brj fix signatures. Special emphasis will 
be given on the LMNDM scenarios. In Section 3 we shall examine the hrj fix signal and 
occasionally the hlj ^t signal in the mSUGRA model using the above data and comment 
on the viability of realizing the LMNDM scenarios in view of the uncertainties in the LHC 
constraints. In Section 4 the analyses of Section 2 and 3 are updated in the light of the 
recent C = lfb~^ data. Our conclusions will be summarized in Section 5. 



2 SUSY signatures and LMNDM in generic models 

The simplest generic model compatible with all LHC data accumulated so far (C = 35 pb^^) 
would be one with all strongly interacting sparticles beyond the reach of 7 TeV experiments 
while all sparticles in the electroweak sector are light. Unfortunately the earlier simulations 
in the context of the 14 TeV run indicate that any signature in the current experiments at 
lower energy is not likely. Thus for a model with non-trivial signatures at this stage of the 
LHC experiment one needs at least one relatively light strongly interacting sparticle. 

In the first model with modest values of tan /3, only the third generation of squarks and 
the sparticles in the EW sector are assumed to be within the reach of the early phases of 
7 TeV run. The hlj ^t signal has already been studied in [28] in such a phenomenological 



^In spirit these models are similar to the simplified phenomenological model considered by the ATLAS 
collaboration with only squarks of the first two generation, the gluinos and the LSP within the reach of the 
7 TeV run (see Fig. 2 in the second paper of [1]). 



scenario . The large mixing in the top squark mass matrix producing a hght mass eigenstate 
provides a quahtative justification. For small or moderate tan/3 the b-squark mass eigenstates 
will be much heavier. Of course a large trilinear soft breaking term At is needed for the above 
mixing. It will also yield rrih compatible with the LEP bound through radiative corrections. 

The current LHC data hardly constrain this model since events from direct low mass 
ti — tl pair production have too little /i!t or mg// [28] to survive the strong cuts currently 
implemented by the ATLAS and CMS experiments for SUSY search. Dedicated searches 
with softer cuts are called for. 

In this model the first two generations of squarks and the gluinos are assumed to be 
beyond the reach of the early stages of the on going LHC experiments due to some yet 
unknown soft breaking mechanism. We further assume for the sake of simplicity that rriQ 
is the common mass of the squarks belonging to the third generation, the sleptons and the 
two neutral Higgs bosons. Similarly mi/2 controls the masses of the electroweak gauginos 
only. If this partially constrained spectrum yields observable signal over a reasonably large 
parameter space it is obvious that more will be available in a totally unconstrained MSSM. 
In fact if the unification of the electroweak gaugino masses at Mq is relaxed, the mass of the 
LSP DM candidate can be even lower [30]. 

We shall fix tan/3 = 10 but take Aq as a variable. Again for simplicity the magnitude of 
Ao is restricted to be less than 1 TeV. While computing the spectrum we have checked that 
no charge-colour breaking minimum of the scalar potential [31] occurs. For each point the 
minimum allowed Aq is determined by the rrih bound from LEP [24]. 

The radiative corrections to the lightest Higgs boson mass {rrih) involve some theoretical 
uncertainties (see, for example, [15, 28] for a brief discussion and references to the original 
works). In view of these uncertainties, if the computed Higgs mass is rrih > 110 ^ for a point in 
the parameter space, that point will be regarded as acceptable although the computed mass 
is somewhat smaller than the direct bound from the direct searches at LEP [24]. Throughout 
this work the pole mass of the top (bottom) quark will be taken as mt{mb) = 173 (4.25) 
[32]. We shall assume that the masses of the lighter chargino, all the sleptons except the tau 
mass eigenstates and the third generation squarks are heavier than 100. This is basically a 
simplified form of the LEP limits. The lighter f mass eigenstate is assumed to be heavier 
than 82 as required by the LEP data. 



^See footnote 4. 



The light sparticle masses and decay branching ratios (BRs)are generated by SUSPECT 
[33] and Sdecay [34] and /i is fixed by the REWSB condition [6]. 

In the parameter space of interest for the on going sparticle searches at the LHC the two 
body decays ti — )■ bxt occurs with almost 100% branching ratio (BR). Moreover in bulk of 
the parameter space the decays xt "^ t^i^? ^t dominate, leading to the brj ^t signature. 
When the latter decays become subdominant or is kinematically suppressed the decay into 
a real W {xt ~^ Xi^) or a virtual W {xt ~^ Xi^* '^ Xiff') open up yielding the hlj ^t 
signature [28]. 

We have generated tit\ pair events at Ecm = 7 TeV using Pythia [35]. The signal brj^T 
has been simulated using the following procedure. Initial and final state radiation, decay, 
hadronization, fragmentation and jet formation are implemented following the standard pro- 
cedures in Pythia. 

In this paper all leading order (LO) signal cross-sections have been computed by CalcHEP 
[36] unless otherwise stated. For any two body final state (except for QCD processes) with 
identical particles or sparticles both the renormalization and the factorization scales are 
taken as, fiR = fip = M, where M is the mass of the particle or sparticle concerned. For two 
unequal masses the scales are taken to be the average of the two. For QCD events the scales 
have been chosen to be equal to vS which is the energy in the parton CM frame, and the 
cross-section is computed by Pythia. All LO cross-sections are computed using CTEQ5L 
parton density functions (PDFs) [37]. 

We have considered the backgrounds from ti, QCD events and W + n-jets events, where 
W decays into all channels, ti events are generated using Pythia and the LO cross-section 
has been taken from CalcHEP which is 85.5 pb. QCD processes are generated by Pythia in 
different pr bin : 25 < pr < 400, 400 < pr < 1000 and 1000 < pr < 2000 , where pT is 
defined in the rest frame of the parton collision. The main contribution comes from the low 
pT bin, which has a cross-section of ~ 7.7-E'+07 pb. However, for other bins (400 < pT < 1000 
and 1000 < pr < 2000), the background events are negligible. 

For W + n-jets events we have generated events with n = 0, 1 and 2 at the parton level 
using ALPGEN (v 2.13) [38]. We have generated these events subjected to the condition 
that Pj. > 20, AR{j,j) > 0.3 and |?7| < 4.5. These partonic events have been fed to Pythia 
for parton showering, hadronization, fragmentation and decays etc. 

The next to leading order (NLO) cross-sections for stop-stop pair production have been 



computed by PROSPINO [39] using the CTEQ5M PDFs. The K-factors are computed by 
comparing with the LO cross-section. The LO cross-sections from PROSPINO agree well 
with CalcHEP for the same choice of the scales. 

The NLO background cross-sections are not known for some backgrounds - in particu- 
lar for the QCD processes. For computing the significance of the signal we conservatively 
multiply the total LO background by an overall factor of 2. 

We have used the toy calorimeter simulation (PYCELL) provided in Pythia with the 
following settings. 

• The calorimeter coverage is |?7| < 4.5. The segmentation is given by At] x A0 = 
0.09 X 0.09 which resembles a generic LHC detector. 

• A cone algorithm with A R= v^Arp~+"A0^ = 0.5 has been used for jet finding. 

• E^r^*j^j^ = 30 and jets are ordered in Et. 

The signal has been selected as follows: 

Lepton Veto: 

Leptons (/ = e, /i) are selected with Pt > 10 and |?7| < 2.4. For lepton-jet isolation we 
require AR{l,j) > 0.5. For the sake of simplicity the detection efficiency of e and /i are 
assumed to be 100%. Events with isolated leptons are rejected. 

b- jet identification: 

We have tagged fe-jets in our analysis by the following procedure. A jet with \r]\ < 2.5 
corresponding to the coverage of tracking detectors matching with a i?-hadron of decay 
length > 0.9 mm has been marked tagged. This criteria ensures that single fc-jet tagging 
efficiency (i.e., the ratio of tagged 6-jets and the number of taggable 6-jets) e^ ~ 0.5 in ti 
events. 

r- jet identification: 

Taus are identified through their hadronic decays producing narrow jets with 1 or 3 
tracks pointing to the jets. We have defined a narrow signal cone of size ARs = 0.1 and an 
isolation cone of size ARj = 0.4 around the calorimetric jet axis. We then require 1 or 3 
charged tracks inside the signal cone with \ritrack\ < 2.5 and P^ > 3 for the hardest track. 
We further require that there are no other charged tracks with P^ > 1 inside the isolation 
cone to ensure tracker isolation. 



The following cuts will be call Set 1. These cuts ensure stop rich signal events while 
rejecting the background efficiently : 

• We have selected events with one tagged h jet [cut 1.1). 

• We have selected events with one tagged r jet {cut 1.2). 

• We have rejected events with isolated lepton [cut 1.3). 

• Events with missing transverse energy {^t) > 70 are selected {cut 1.4)- 

• We have also demanded events with P^ tagged r jet > 40 [cut 1.5). 

or 

We have demanded events with Px tagged b iet < 50 {cut 1'.5). This cut is particularly 
useful if the mass difference between the ti and the xf is small. It also rejects the ti 
background efficiently. 

The set of cuts which includes cut 1' .5 will be referred to as Set 1' . 

Our results are summarized in Fig. 1. In the green (online) or small crosshatched (red 
(online) or big crosshatched)) region the hrj fix signal can be observed for an integrated 
luminosity of £ = 1 (5) fb~^. For observability we simply require S/\fB > 5, where S (B) is 
the number of signal (background) events. We find that mj upto 235 (280) can be probed 
by this signal with £ = 1 (5) fb~ . 

In the blue (online) or hatched region the decays xt ~^ ^i^n ^r^ are either phase space 
suppressed or are kinematically forbidden. But the hlj ^t signal becomes observable with 
the cuts of Set 1 proposed in [28] for a suitable Aq. In the blue (online) or hatched region 
mj varies between 130 - 190 for suitable choices oi Aq. Beyond a certain ttiq (approximately 
480), Aq > 1 TeV will be required to make mj sufficiently light for an observable signal. 

To get a feeling for the relative size of the signal and the SM backgrounds we consider 
four benchmark points BPl (110,170,-540,10), BP2 (110,190,-640,10), BP3 (130,200,-550,10), 
BP4 (150,230,-745,10) which yield observable signals at £ = 1 fb"^ (see Fig. 1). The 
quantities in the bracket are rrao,mi/2, ^o a-nd tan/3. The stop mass at the four points are 
143.6, 165.0, 215.6 and 208.2 respectively. At all four benchmark points BR(ti — )> hxt) and 
BR(xi^ -^ Tiv) are maximal to a very good approximation. 
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The 0/ + jets + ^t or any other signal which gets contribution from all squark-gluino 
production is not particularly sensitive to the presence or absence of a light stop. However 
the brj^T signal proposed here stems from ti—tl pair production alone. Hence it sensitively 
depends on mj and , consequently, on Aq. The quoted Aq gives the largest signal in each 
case. However, for each point a range of Aq values leads to observable signals. In this range 
the ti-Xi mass difference corresponding to |Ao|max is too small for producing a taggable b-jet 
and the signal becomes weak. In contrast at |Ao|mm, the ti — tl production is too suppressed 
- due to a relatively heavy ti- to produce a viable signal. This is true for all points in Fig. 
1. For example with rriQ = 90, mi/2 = 180, the signal is observed in the range l^ol = 370 - 
550. For e.g., Aq = -370, -500 and -550 observable signals occur for £ = 5 fb~^, £ = 1 fb~^ 
and £ = 5 fb~^ respectively. 

The LO cross-sections after different cuts for the four benchmark points and different 
SM backgrounds are presented in Table 1. 

The last row gives the significance of the signal for £ = 1 fb~ for the two sets of 
cuts. The significances, however, are computed on the basis of the NLO cross-sections as 
discussed above. For the range of m-jr relevant here the K-factor varies rather slowly. We 
have multiplied the LO signal cross-sections by the average value which is approximately 
1.7. As stated above the total SM background is multiplied by an overall factor of 2. 

It follows that larger significance is obtained in some regions by using cut 1.5 while in 
others cut 1' .5 is more appropriate. The cut 1.5 will be more effective in the region where 
the P^ is harder. This is the case for BP3 as a result of the relatively large mass difference 
between fi and Xi- On the other hand for regions where the difference between mj and 'rn~+ 
is less than 50, cut 1' .5 will be more effective which is the case for BPl, BP2 and BP4. 

In Fig. 1 each point marked with a cross is consistent with the DM relic density data ^. 
In Figure 1 we also identify the parameter space sensitive to the proposed signal. Therefore, 
in this reach plot we do not restrict ourselves to points allowed by the WMAP data only. 
In fact at many points in the delineated parameter space the computed Vth"^ violates even 
the WMAP upper bound on the relic density. If a signal is indeed observed at any of these 
points, it would indicate that the lightest neutralino though stable in the time scale of a 



^We remind the reader that the parameters toq and mi used here are different from the common scalar 
and gaugino masses in mSUGRA. We have defined these parameters in the present context in Section 2, 
paragraph 4. 
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Signal 


Background 


a(pb) 


BPl 

38.9 


BP2 

19.2 


BP3 

4.56 


BP4 

5.54 


ti 
85.5 


QCD 

7.7 X 10^ 


W + lj 
1.43 X 10^ 


W + 2j 
5200 


cut 1.1 


4.6251 


3.5190 


2.2193 


1.6781 


42.1780 


6.454 X 10^ 


44.0797 


48.1231 


cut 1.2 


1.0696 


0.6797 


0.7327 


0.5540 


3.8987 


6.426 X 10^ 


2.7408 


2.9982 


cut 1.3 


0.8581 


0.5443 


0.5365 


0.4164 


2.8266 


6.348 X 10^ 


2.4310 


2.6816 


cut 1.4 


0.2899 


0.1887 


0.2469 


0.1689 


0.7301 


22.1066 


0.1787 


0.3067 


cut 1.5 
or 

cut 1'.5 


0.2075 
0.2067 


0.1262 
0.1273 


0.1802 
0.0784 


0.1278 
0.1045 


0.5318 
0.1855 


0.0102 
0.0001 


0.1191 
0.0595 


0.2126 
0.0940 


s/Vb 

After cut 1.5 
After cut 1' .5 


8.4 
13.4 


5.1 
8.3 


7.3 

5.1 


5.2 
6.8 











Table 1: The LO cross-sections (including efficiency) for the signal corresponding to BPl - 
BP4 in the phenomenological model ( section 2, paragraphs 2-4) and the SM backgrounds 
after Set 1 or Set 1' of cuts. The last row gives the significance of the signal for £ = 1 fb"'^ 
on the basis of NLO cross-sections (see text). 

collider event may be cosmologically unstable. A tiny R-parity violation induced by higher 
dimensional operators, for example, may induce neutralino decays at the cosmological scale. 
In such cases the observed relic density must come from some other non-neutralino sources. 

In [15] it was shown that there are two generic regions corresponding to LMNDM. There 
is a region where neutralino pair annihilation via R-type light slepton exchange or bulk 
annihilation produces the observed relic density. The tension between the computed lighter 
Higgs scalar mass (m/j)and the corresponding experimental lower bound is softened by non- 
zero but moderate negative values of Aq (200 - 300) and the uncertainties in the computed 
771/j. Here the brj^T signal is observable since the fi - Xi mass difference is relatively large. 
Hence taggable r-jets arise from the decay fi — t- tx^. 

In the other regions consistent with LMNDM the lighter f mass eigenstate is much lighter 
due to larger values of \Aq\ [O (1 TEV)). As a result fi - xH coannihilation along with bulk 
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Figure 1: Regions of mo - ''tii/2 plane which can be probed by direct production of titl 
pairs for 1 < £ < 5 fb~^ (see Section 2, paragraph 2 -4 for the details of the underlying 
phenomenological model). The SM backgrounds are in Table 1. 

annihialion produce the LMNDM. However, due to the small fi - LSP mass difference the 
brj^T signal from stop pair production may not be viable. 

The next scenario we consider is one with two strongly interacting sparticles within 
the reach of the first phase of the LHC experiment along with electroweak sparticles. It is 
assumed that the production of the lighter top squark and the gluinos is the primary source 
of the LHC signatures. The other strongly interacting sparticles are assumed to be beyond 
the reach of the LHC-7 TeV run. The masses of all strongly interacting sparticles are taken 
to be independent of the electroweak sector. For simplicity the masses of the sparticles and 
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other relevant parameters in the EW sector are chosen as in mSUGRA. One is thus free 
to choose any LMNDM scenario. This will be referred to as the Light Stop Gluino (LSG) 
scenario. 

We shall first consider a phenomenological model with ?n^ and m^ as unrelated param- 
eters at the weak scale. We shall also consider a mSUGRA type scenario with non-universal 
boundary conditions at the GUT scale which leads to the LSG scenario. 

However, the parameter space in the LSG scenario is already constrained to some extent 
by the 35 pb~^ LHC data. The ATLAS collaboration has analysed the jets + JUt signal 
using four sets of cuts referred to as A, B, C, D [1]. The corresponding lower limits on the 
production cross-sections of all strongly interacting sparticles including the efficiencies are 
1.3, 0.35, 1.1 and 0.11 pb respectively. These limits are converted into constraints in the 
777,0 ~ ^^1/2 plane by computing the NLO cross-section by PROSPINO and the efficiencies by 
Pythia. 

In mSUGRA the set D consisting of very hard cuts is very potent for obtaining new mass 
bounds. In the LSG scenario, however, this set kills the signal in any LSG model consistent 
with the Tevatron bounds on m^ for heavy squarks [4]. When each point in the parameter 
space is required to pass all the above cut sets, set A and set C turns out to be most effective 
in obtaining limits in the LSG scenario. The resulting bounds are m-x > 320-330 for mj = 
150 - 300. The insensitivity of the bound on m^ to mj clearly indicates that the hard cuts 
employed by the ATLAS group in isolating the jets + ^t signal eliminates all events from 
stop pair production. The bound essentially comes from gluino pair production followed by 
the decay of each gluino into top-stop pairs. Hence the canonical signals with strong cuts 
are not suitable for revealing the presence of a light stop. Interesting contribution to the 
signal may, however, originate from gluino decays. 

For computing the above limits we have set the masses of the squarks (also slep- 
tons) of L and R type belonging to the first two generations at 1.5 TeV. We have treated 
TTij- , mj , m-j^ ,At,Ai) as well as m^ as free parameters. We have chosen the following param- 
eters for the electroweak sector at the weak scale as Mi =60, M2 =125, m~^ =155, m~^ =116, 
At- = —615, tan/3 = 10, /i = 348 . This choice of parameters gives m~± = 122, m~o = 59, 
?7i~ = 109, m~^ = 142 and compatible with the WMAP data. The above limits mildly 
depend on the variation of the parameters in the EW sector unless one makes very specific 
choices such that the decay pattern of the EW gauginos are drastically different ( e.g., they 
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all decay leptonically with 100% BR). 

SUSY searches at the Tevatron can also potentially constrain the LSG scenario. We have 
not done the analysis in this paper. As a cautious approach we have only considered m^ ~ 
400 which corresponds to the bound on m^ if all squarks are heavy and is stronger than the 
limits discussed above. 

We next compute the hrj ^t signal for different mj and m-x in the LSG scenario. The 

1.5 TeV, m~± = 152.6, m~o 



other parameters are fixed as follows: rra, 



gL/«(i.2) 



m- 



II 



m-r 



= 79.1, rra~ = 144.1. We take wino like Xi ^ind bino like Xv The electroweak sector is 
chosen such that the resulting DM relic density is consistent with WMAP data. 



Points 


htl 


LSG scenario 


mj^-m~ 


Cut Set 1 


Cut Set 2 


Cut Set 3 


Cut Set 1 


Cut Set 2 


Cut Set 3 


209 - 504 


359.0(8.6) 


270.3(8.8) 


3.8(0.9) 


463.4(11.1) 


280.5(9.1) 


28.7(7.2) 


250 - 443 


204(4.9**) 


136.7(4.4**) 


6.8(1.7) 


285.3(6.8) 


144.7(4.7**) 


25.4(6.3) 


250 - 517 


209(5.0) 


140.7(4.5**) 


6.6(1.6) 


239.2(5.7) 


139.2(4.5**) 


20.1(5.0) 



Table 2: Number of hrj ^t events from pure tit\ production (columns 2-4) and that in the 
phenomenological LSG scenario (columns 5-7) for three sets of cuts (see text). Numbers 
in the brackets are the significance of the signal for C = 1 fb~^ using NLO cross-section and 
the entries marked with ** indicate that the signal is observable for 1 < £ < 5 fb~^. 

We analyse the signal with 3 sets of cuts and the results are summarized in table 2. The 
Cut Set 1 and the corresponding SM background have already been given. Both stop and 
gluino production can contribute to the signal. If the stop is relatively heavy and the gluino 
is relatively light the size of the overall BSM signal may improve considerably to the case 
where only ti of the same mass is present. 

The Cut Set 2 has two cuts in addition to Cut Set 1: 

• Ncentrai-jet < 4 ( where the central jets have \ri\ < 2.5) . 

• Meff < 500 ( where Me// = \f}r\ + ^i\P^'\ + ^ilPrl {h = e,/i)). 

The corresponding SM background is 0.4744 pb. After this set of cuts the gluino induced 
events in the sample are drastically reduced. As discussed in [28] this sample can be used 
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for studying the properties of the hght stop squarks. Moreover, if the signal conies from stop 
pair production alone the fraction of events surviving this cut is much larger than that if the 
signal stems from both ti and gluino induced events. This distinction , however, becomes 
unclear as m^ increases. 

The Cut Set 3 is designed to remove the pure stop induced events. It consists of the 
following stronger cuts : 

• We have selected events with one tagged b jet {cut 3.1). 

• We have selected events with one tagged r jet {cut 3.2). 

• We have rejected events with isolated leptons {cut 3.3). 

• Events with at least 6 central jets are selected, where central jets are defined as pycell 
jets with \ri\ < 2.5 {cut 3.4). 

• Events with missing transverse energy {^t) > 160 are selected {cut 3.5). 





ti 


QCD 


W + lj 


W + 2j 


cut 3.1 


42.1780 


6.454 X 10^ 


44.0797 


48.1231 


cut 3.2 


3.8987 


6.426 X 10^ 


2.7408 


2.9982 


cut 3.3 


2.8266 


6.348 X 10^ 


2.4310 


2.6816 


cut 3.4 


0.4454 


0.0256 


0.0476 


0.2275 


cut 3.5 


0.0076 


0.0004 


- 


- 



Table 3: The LO cross-section (including efficiency) of SM backgrounds after the Set 3 of 
cuts. 

Only The gluino induced events survive and depending on the gluino mass may give an 
observable signal. The significance of the signal for £ = 1 fb~^ for each case is given in 
parentheses in Table 2. The corresponding SM backgrounds are presented in Table 3. For 
the gluino mass range in Table 2 the K factor varies between 1.6 and 1.7. 

The LSG scenario can be realized in gravity mediated SUSY breaking with non-universal 
masses at Mq. All scalar superpartners squarks belonging to the first two generations are 
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assumed to have masses beyond the reach of the ongoing runs at the LHC The squarks 
belonging to the third generation are assumed to have much smaller masses. At Mq the 
T mass is assumed to be even smaller so that the LMNDM can be eventually realized. 
Qualitatively the non-universal masses for the third generation can be generated by the 
running of the corresponding soft breaking parameters, from the SUSY breaking scale down 
to Mg [40], by effective mass terms induced at Mq by some yet unknown flavour dependent 
interactions above the GUT scale etc, etc. The phenomenology of these models has attained 
due attention [41]. 

The gluino mass (M3) is assumed to be smaller than the other gaugino masses at Mq- 
For simplicity it is assumed that Mi = M2 at this scale. They are chosen such that they are 
compatible with a LMNDM scenario. The motivation for non-universal gaugino masses have 
already been discussed by several authors [42]. It arises if a GUT nonsinglet chiral superfield 
couples to the gauge kinetic function and the hierarchy among the non-universal gaugino 
masses depends on the representations of the GUT group to which the chiral superfield 
belongs. In practice linear combinations of such chiral superfields may couple, making the 
prediction of the above mass hierarchy rather difficult. The phenomenology non- universal 
gaugino masses has also been discussed extensively [43]. Models with both the above non- 
universalities have also been constructed [44] and the resulting phenomenology were analyzed 
[45]. 

Finally in order to get the magnitude of /i consistent with a LMNDM scenario with bino 
like LSP, we have to introduce non-universal soft breaking Higgs masses ttih^ and trh^ at 
Mg [46]. 

For illustration we have chosen the following spectrum at Mg which is consistent with 
a LMNDM scenario satisfying the WMAP constraints, where the relic density is produced 
by stau - LSP coannihilation: mo(l,2) = 1.5 TeV, mH^= 300, mH^= 500, Mi = M2= 240, 
At = Ab = Ar = -500, tan/3 = 10. 

The parameter mo(3)(t), the common mass for the third generation of squarks {Ql-, tn, 
6ij) and the gluino mass M3 at Mg are the variables. The common mass ?Tio(3)(r) of the 
third generation of sleptons (r^, r^, z/^) is mildly varied in a small range to obtain a LMNDM 
scenario. The variation of the soft breaking terms at Mg and the resulting physical masses 
of the relevant sparticles are presented in Table 4. 

In the phenomenological LSG scenario mj is a weak scale input and At is not important. 
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Ms 


mo(3) 
it) 


mo (3) 

(r) 


^t. 




m~ 


m~o 

Xi 


rUg 


SI 
{Cut Set 1) 


S3 

{Cutsets) 


150 


460 
520 
580 


112 
117 
121 


214 
288 
351 


187 
190 
192 


106.3 

107.4 
107.7 


98.0 
98.6 
98.9 


429 
435 
438 


329.8(7.9) 
196.3(4.7**) 
136.1(3.3**) 


35.2(8.8) 

20.1(5.0) 

18.7(4.7**) 


170 


460 


114 


250 


188 


106.6 


98.0 


475 


235.6(5.6) 


23.6(5.9) 


200 


460 


117 


302 


189 


107.4 


98.0 


547 


75.8(1.8) 


10.6(2.6**) 



Table 4: Some representative non-universal input parameters at Mq leading to the LSG 
scenario and the corresponding sparticle spectra. The EW sector yields DM relic density 
consistent with WMAP data. Last two columns give number of hrj ^t events (using NLO 
cross-sections) under cut Set 1 and Set 5 for £ = 1 fb~ . Significance for each case is given 
in the bracket and entries marked with ** indicate that signal is observable for 1 < £ < 5 
fb-\ 

In the LSG scenario obtained from non-universal GUT scale boundary conditions, a smaller 
mj contributes more to the size of the hrj J^t signal obtained with softer cuts (see the last 
but one columns of table 4 and section 4). Hence At at Mq is an important parameter (see 
the inputs for table 4). 

The prospect of observing the 6rj fjT signal at the ongoing LHC experiments in the above 
LSG scenario is also summarized in Table 4. It follows that for a fixed M3 = 150 (m^ ~ 435), 
vTLj < 350 can be probed with £ a little more than 1 fb^^. 

The hrj ^^ signal, which arises when the chargino decays into r rich final states with 
large BR, has not so far been searched by the LHC collaborations. The ATLAS collaboration 
has already searched for the hlj J^t signal arising from the decay of stop pairs using the 35 
pb~^ data (see the third paper in ref. 1). Such pairs may either be produced directly or via 
two body decay of gluino pairs produced at the LHC. Apparently the experiment takes into 
account both possibilities and excludes m-x below 530 - 540 for mq = 125 - 300. However, the 
insensitivity of the gluino mass limit on mj indicates that the events from directly produced 
stop pairs are eliminated by the selection procedure and effectively their limits arise from the 
gluino induced processes. As already shown in [28] and also in this paper directly produced 
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light top squarks can also contribute to the signal under a different search strategy based on 
softer cuts. The possibility that only lighter stop squark pairs have been produced at the 
LHC yielding either the blj ^t or brj ^t signature still remains open. 

We first employ relatively soft cuts so that the search is sensitive to both direct stop pair 
production and gluino pair production followed by the decay of each gluino into top-stop 
pairs. This enhances the size of the overall SUSY signal especially for relatively heavy stop 
{mj > 250) and light gluino. As demonstrated above one can design cuts which can separate 
the signal stemming from the two different production channels. 

As already noted the simplest model practically unconstrained by the LHC data would 
of course be the one in which all strongly interacting sparticles are beyond the reach of 
the LHC experiments at 7 TeV. This model can be motivated by introducing non-universal 
gaugino masses in a mSUGRA like framework by requiring M3 >> mi/2, rriQ at Mq , where 
M3 is the gluino mass and mi/2 is the common mass of the SU{2) and U{1) gauginos 
[47]. The rest of the sparticle masses may be determined by mSUGRA like parameters 
mo, Ml = M2 = ?7ii/2,^0; tan/3 and sign (/i). The parameter mo, assumed to be much 
smaller than M3 at Mq, controls the scalar masses in the electroweak sector but has very 
little impact on the masses of the strongly interacting sparticles. The latter parameters can 
be chosen such that they are consistent with the LMNDM scenario. 

Unfortunately the signal will not be easy to observe at the on going LHC experiments. 
The slepton pair production [48] and the clean trilepton signal [49] from chargino (xi')-second 
neutralino (X2) pair production would be the most distinctive signatures of this model. It is, 
however, well known that even for the 14 TeV experiment the reach in the trilepton channel 
is rather modest [50]. This does not suggest an exciting prospect for the 7 TeV run. At 14 
TeV the reach can be improved by including the 1/ + 2 r or 2/ + Ir events arising from the 
decays of xt ^.nd X2 [47]. The reach for slepton search is also rather limited at 14 TeV [51] 

These signals, however, have not been studied for the on going run nor is the prospect 
of slepton search at 7 TeV known. In this paper we did not pursue this model any further. 
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3 SUSY signatures and LMNDM in mSUGRA in view 
of LHC data with £ = 35 pb"^ 

In this section we consider the prospect of observing the brj^x signal at LHC-7 TeV and re- 
ahzing the LMNDM scenario in mSUGRA. However, the ATLAS collaboration has obtained 
strong constraints on the mSUGRA parameter space [1] which should be taken into account. 
We remind the reader that the ATLAS results were obtained for Aq = and tan/3 = 3 - a 
choice forbidden by the direct bound on rrih from LEP. However, in the same paper they 
noted that the limits are rather insensitive to the variation of the above two parameters. We 
shall, therefore, assume that these limits are valid in more general models with different Aq 
and tan/3 values. 

Our results are summarized in Fig. 2. This figure corresponds to Aq = —600 and tan/3 = 
10. The observed ATLAS limits (reproduced in Fig. 2) are rather strong and strongly 
disfavors the LMNDM scenario. However, these limits were obtained with the strongest set 
of cuts used by them (signal region D in [1]) resulting in 2 observed events against a SM 
expectation of 2.5 ± 1.0^o4 ± 0.2. An upward fluctuation of the observed number, which 
is certainly a distinct possibility, would relax the limits significantly. Thus the variation 
of the expected limit, which is very close to the observed limit in this case, in response to 
±la fluctuation of the SM expectation is perhaps a more realistic description of our present 
knowledge. The band resulting from the above fluctuations as given in [1] is also included 
in Fig. 2. 

We also present the points (denoted by the x mark) in the low mQ-mi/2 region allowed by 
the DM relic density data. The relic density is mainly produced by stau-LSP coannihilation 
with some contribution from bulk annihilation. The observed limits disfavor most of the 
region consistent with the LMNDM scenario. In fact m~o < 143 is disfavored. However, if 
the uncertainties in the limits are taken into account Xi with smaller masses can not be the 
excluded with certainty as DM candidates. It is also to be noted that the computed rrih for 
Aq = also puts pressure on the low mo - ^^1/2 region even after theoretical uncertainties in 
the computation is taken into account. 

The Cut Set 3 introduced in Section 2 is used for estimating the signal and the back- 
ground (see Table 3). In the green (online) or small crosshatched region the signal is observ- 
able with £ = 1 fb~ . No signal from the points consistent with the observed relic density is 



20 



550 



500 



450 



400 



350 



300 



250 



200 



150 




100 



200 



300 



400 



500 600 

rrin 



700 800 900 1000 



Figure 2: The green (online) or small crosshatched (red (online) or big crosshatched)) region 
of 777,0 - ^^1/2 plane in mSUGRA can be probed by the brj ^t signal from all squark-gluino 
events using Set 3 of cuts with £ = 1 (5) fb~^. The blj JUt signal probes the grey (blue 
(online) or hatched) region with the cuts proposed in [28] (cut Set 4 introduced in this 
paper) for 1< £ < 5 fb~^. Here Aq = -600 and tan/3 = 10, sign(/i) > (see text for the 
details). 

expected for DM allowed points even after considering the uncertainties in the ATLAS data. 
For 1 fb~^ < £ < 5 fb~^ the signal is observable over a much larger region of the parameter 
space (the red (online) or big crosshatched region). A few points allowed by the relic density 
data yield observable signals (see fig. 2). 

We next comment briefly on the blj J^t signal. In addition to the Set 3 of cuts defined 
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in [28] we introduce a new set of cuts [Set 4) given below. 

• We have selected events with one tagged b jet. 

• We have selected events with one isolated lepton . 

• Events with at least 6 central jets are selected, where central jets are defined as pycell 
jets with |?7| < 2.5. 

• Events with missing transverse energy (^t) > 160 are selected. 

The former set gives a better reach in mi/2 for rriQ < 550 while the second set is more effective 
for larger ttt-q for £ = 5 fb^ . This is illustrated in Fig. 2 by the grey and the blue (online) 
or hatched regions respectively. 

Recently several groups have reexamined the sensitivity of the on going XENONIOO 
experiment in the light of the constraints on the mSUGRA model imposed by the LHC 
35 pb" data. The Figure 2(L) of [20], for example, indicates that the parameter space 
sensitive to the above experiment is not compatible with the LHC constraints. However, if 
the sensitivity of XENON experiment is improved by a factor of 10, which may be achieved 
by the end of 2012, a much larger region of the mSUGRA parameter space may be accessible 
to direct search experiments. Discovery of neutralino DM by direct search in the near future 
would, therefore, suggest models beyond mSUGRA like the ones discussed in this paper. 

4 The generic models, the mSUGRA model and the 
LMNDM scenario revisited in view of the C = 1 fb~^ 
data 

After submitting this paper to the arXiv and the journal, the constraints on the mSUGRA 
model from LHC data with £ = 1 fb~ were announced. In this section we discuss briefiy 
the impact of the new data on our main results. The following discussion is based on the 
results presented in the Lepton-Photon conference, 2011 [52] and in the CMS paper [53]. 

Of course the results for the brj /^j- signal from stop pair production alone in the first 
part of Section 2 remain unaltered. As already discussed, none of the CMS and ATLAS 
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searches employing very hard J^t cuts are sensitive to either the blj /^t signal discussed in 
[28] ^ or the brj /^t channel discussed in this paper. A dedicated search for relatively low 
mass top squarks using softer cuts as out lined in Section 2 is called for. Also the search for 
final states with one tagged b-jet and one r jet has not been reported so far by the LHC 
collaborations. 

The stronger constraints from the new data in the jets + ^t channel will increase the 
lower bounds on m-^ in the LSG scenario presented in Section 2. However, the details of the 
cross section limits similar to those presented by the ATLAS group for the 35 fb~^ data used 
in obtaining our bounds, are not yet available for the £ = 1 fb~^ data. We, therefore, use 
the following procedure. 

From the ATLAS analysis in the jets + ^t channel based on 165 pb^^ of 2011 data [55] 
the limit on gluino mass is 500 for heavy squarks(?Tig = 1250). The CMS analysis in the 
same channel using 1.14 fb~^ of 2011 data [53] puts a stronger limit on the gluino mass. For 
heavy squarks with average mass 1.5 TeV (2 TeV), the lower limit on gluino mass is 570 
(550). In the Lepton-Photon conference, 2011 [52] ATLAS results in the same channel was 
presented for 1.04 fb~^ of 2011 data. For heavy squarks gluino mass below 600 are excluded. 

Strictly speaking the above limits are not directly applicable in the LSG scenario. In 
this case a gluino decays via the cascade ^ — ?■ tit — ?■ {bx^){bW^), yielding the signal whereas 
CMS and ATLAS consider gluino pair production in the limit where all squarks are heavy, 
so that each gluino decays via 3 body modes(ggx'j') QQX2 or QQ'x"^)- Jets are, therefore, likely 
to be softer in our case on the average. Hence the hard cuts on jet p^ employed by the LHC 
collaborations are expected to give somewhat weaker limit on m-^ in the LSG scenario. As a 
reasonable guess we have considered m^ > 550 in the LSG scenario. Results for m^ < 500 
are already available in Section 2. 

The ATLAS group has also updated their search for the blj /^t signal [54] in the LSG 
scenario. Using 1.03 fb~ data they have excluded gluino masses below 500-520 mj in the 
range 125 - 300. This limit, however, is also not applicable to the analysis in this paper 
based on a different chargino decay mode (BR(x^ — ?■ fiZ/T-) = 100 %). 

We present in Table 5 the observability of the signal in the phenomenological LSG model 
where m-^ and m-x are unrelated parameter. This table illustrates the search prospect for 



^The above comment is also applicable to the negative search in blj^T channel by the ATLAS collabo- 
ration [54] briefly reviewed below. 
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Points 


titl 


LSG scenario 


^t. 


m~ 


Cut Set 1 


Cut Set 2 


Cut Set 3 


Cut Set 1 


Cut Set 2 


Cut Set 3 


209 


550 
700 
850 


362.9(8.7) 


272.3(8.9) 


4.5(1.1) 


394.4(9.4) 
378.3(9.0) 
374.2(8.9) 


285.8(9.3) 
275.4(9.0) 
270.3(8.8) 


17.4(4.3**) 
10.3(2.5**) 

8.8(2.2)** 


250 


550 
650 

750 


207.6(5.0) 


136.0(4.4) 


6.8(1.6) 


230.7(5.4) 

210.8(5.1) 

204.3(4.9**) 


138.2(4.4**) 
136.5(4.4**) 
132.1(4.2**) 


16.3(4.0**) 

10.2(2.5**) 

8.3(2.1) 



Table 5: Notations, conventions and input parameters are the same as Table 2. 
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mo(3) 
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mo (3) 
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"^Ii 
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m~ 


m~o 


m~g 


SI 
{Cut Set 1) 


S3 

{Cut Set 3) 


210 


365 

420 


112 
114 


229 
283 


186 

188 


106 
105 


97 
98 


567 
569 


228.7(5.5) 
101.9(2.4**) 


16.7(4.1**) 
12.2(3.0**) 


225 


330 


112 


231 


185 


106 


97 


600 


218.8(5.2) 


13.9(3.5**) 


250 


250 


111 


235 


184 


106 


97 


657 


201.3(4.8**) 


9.8(2.4**) 



Table 6: Notations, conventions and input parameters are the same as in Table 4. 

'^o ^ 550. The notations, input parameters consistent with a LMNDM scenario and cuts 
are the same as in Table 2. The corresponding results for the scenario with nonuniversal 
boundary conditions at the GUT scale (see Section 2) are presented in Table 6. It may also 
be noted that if £ = 10 fb^^ is accumulated in the LHC 7 TeV run then mass reach in this 



channel will be considerably improved. For example if m 



9 



550 then mr ~ 
ti 



315 can be 



probed with Cut Set 1. 

The new data from LHC based on £ = 1 fb"'^, however, strongly disfavor the LMNDM 
scenario in the mSUGRA model irrespective of the experimental uncertainties. We super- 
impose on Figure 2 the region excluded by the CMS collaboration [53]. It is readily seen 
that for low mo, neutralino - stau coannihilation is the only viable mechanism for DM relic 
density production provided m~o > 215. 
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Even if the issue of neutralino dark matter is set aside, the brj ^t is unobservable in the 
mSUGRA model even for £ = 5 fb~ (see the recent CMS constraint reproduced in figure 
2). The blj /^T signal is observable in a small corner of the parameter space thanks to the 
experimental uncertainties as reflected by the CMS expected limit ±lcr band. 

5 Summary and conclusions 

The data from the on going experiments at the LHC has put a question mark on the via- 
bility of a substantial region of the mSUGRA parameter space corresponding to low mass 
sparticles. This leads to a tension between the data and the low mass neutralino dark matter 
(LMNDM) scenario where neutralino annihilation and/ or neutralino - slepton coannihilation 
can produce DM relic densities consistent with the WMAP data. Moreover, large regions of 
the parameter space sensitive to the ongoing experiments for direct DM search or searches 
in the near future are also under pressure. 

However, the LHC experiments are primarily sensitive to the masses of the strongly 
interacting sparticles - the squarks and the gluinos. In contrast in typical LMNDM scenarios 
the relic density may depend entirely on the properties of the sparticles in the electroweak 
sector. Thus the above tension is an artifact of the model dependent correlations among 
the soft breaking masses in the strong and the electroweak sectors of mSUGRA. Any model 
in which the masses in the two sectors are independent parameters could be free from this 
tension. 

In view of this we propose a few generic models which are unconstrained or mildly 
constrained by the LHC data. These models are generic in the sense that their acceptability 
depends on certain mass hierarchies in the strong sector and not on specific mass values. 
The electroweak sector is assumed to be independent of the strong sector. In fact we only 
assume that the LSP is bino like, tan /3 has intermediate values and all parameters in the 
EW sector are consistent with the corresponding LEP limits provided such limits are not 
based on mSUGRA dependent assumptions. 

In the first phenomenological model under consideration all strongly interacting sparticles 
except for the lighter top squark are assumed to be beyond the reach of the experiments 
at LHC-7 TeV (see Section 2). This squark and the electroweak sparticles are the only 
sources of SUSY signals. For simplicity we assume that the masses of the relatively light 
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sparticles are correlated as in mSUGRA. Our conclusions will obviously hold in a more 
general framework. In this case the trilinear soft breaking parameter [Aq) must be non-zero. 
This ensures consistency with the bound on the Higgs mass from LEP [15]. Moreover, for a 
given rriQ, LSP-fi coannihilation may produce the observed relic density for values of 771,1/2 
significantly lower than that for Aq = [15]. This also facilitates the LMNDM scenario. 

The conventional SUSY signals with hard cuts on /Ey or rrieff are not viable in this 
model. However, a fairly large region of the parameter space can be probed by the hrj J^t 
signal (see Table 1 and the reach plot Fig. 1) using the search strategy sketched in this paper. 
It is estimated that mj < 235 (280) can be probed with £ = 1 (5) fb~^. A part of this 
parameter space is consistent with the WMAP data where the DM relic density is produced 
by bulk annihilation as shown in Fig 1. However, this result is subject to certain simplifying 
assumptions (see Section 2 paragraph 4). In an unconstrained MSSM the signal will be 
compatible with a much larger parameter space consistent with the LMNDM scenario. 

The complementary signal blj ^t signal proposed in [28] may be useful if the charginos 
do not dominantly decay into tau rich final states. The possibility that signals from stop 
pair production only are already buried in the LHC data is still open. 

We next focus on a model (the LSG scenario)in which only ti and g but no other 
strongly interacting sparticle, are within the striking range of the ongoing LHC experiments. 
As before the EW sector is assumed to be independent of the strong sector and is taken 
to be consistent with the WMAP data. This model is only mildly constrained by the LHC 
data for £ = 35 pb"^(see Table 2 for some representative examples). The brj ^^ signal can 
probe a fairly large region of the parameter space (Table 2). Moreover, using additional cuts 
the stop induced events can be separated from the gluino induced ones (see Table 2). The 
sample thus separated can be used for reconstruction of the properties of ti. 

The very recent data for C = lfb~^, however, make the lower limit on m^ more stringent 
(rrig > 550) for heavy squarks (see Section 4). Thus the contribution of gluino pair produc- 
tion to the signal reduce significantly. Still the stop induced and gluino induced contributions 
can be separated by suitable cuts as illustrated in (Table 5). 

The LSG scenario can be realized in a mSUGRA type model with nonuniversal boundary 
conditions at Mq (see Table 4 for sample results for £ = 35 pb~ ).The electroweak sector 
chosen is consistent with a LMNDM scenario. The updated results corresponding to £ = 1 
fb~^ data are presented in Table 6. 
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Finally we look into the brj^T signal in the mSUGRA model (Section 3, Figure 2). We 
also revisit the hi j fir signal in case it gives a better reach in some parameter space. So 
far as the £ = 35 pb~ data is concerned both the above signals and the LMNDM scenario 
cannot be strictly ruled out due to the uncertainties in the data with low statistics ( the 
ATLAS constraints are reproduced in Fig. 2 for ready reference). However, the more recent 
data for C = lfb~ strongly disfavors both the signals and the LMNDM scenario (see the 
recent CMS constraints superimposed on Fig 2). 

The latest LHC data also disfavors a large parameter space in mSUGRA with low mass 
neutralinos in mSUGRA sensitive to the direct DM search by the XENONIOO experiment. 
Alternative models for LMNDM may, therefore, call for more attention. 

In this paper we have not considered constraints from flavour physics like the flavour 
violating decays of the B-hadrons. Strictly speaking these constraints are sensitive to the 
additional assumption that the quark and the squark mass matrices are aligned in the flavour 
space so that the same matrix as the CKM matrix also operate in the squark sector. This 
assumption of minimum flavour violation fails even if there are small off-diagonal elements 
of the squark mass matrix at the GUT scale. On the other hand such small element does 
not affect processes like neutralino annihilation and squark-gluino production and decay. In 
fact, it has been explicitly shown that such small mixings at the GUT scale can significantly 
weaken the constraints from flavour physics [56]. Moreover a comprehensive analysis of all 
possible constraints on a SUSY model is beyond the scope of this paper. 

We conclude that any model whose strong sector is not sensitive to the current LHC data, 
may be consistent with a LMNDM scenario satisfying the WMAP constraints thanks to the 
properties of the EW sector, provided model dependent correlations among the parameters 
of the strong and the electroweak sectors, as in mSUGRA, are given up. Thus there is no 
serious conflict between the LMNDM scenario and the LHC data. 
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